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ABSTRACT 

From 2001 January to 2002 June, we monitored PSRs B0329+54, B0823+26, 
B1929+10, B2020+28 and B2021+51 using the Nanshan 25-m radio telescope of 
Urumqi Observatory to study their diffractive interstellar scintillation (DISS). The 
■ average interval between observations was about 9 days and the observation duration 

ranged between 2 and 6 hours depending on the pulsar. Wide variations in the DISS 
, parameters were observed over the 18-month data span. Despite this, the average 

scintillation velocities are in excellent agreement with the proper motion velocities. 
The average two-dimensional autocorrelation function for PSR B0329+54 is well de- 
. scribed by a thin-screen Kolmogorov model, at least along the time and frequency 

axes. Observed modulation indices for the DISS time scale and bandwidth and the 
^ <— | pulsar flux density are greater than values predicted for a Kolmogorov spectrum of 

O electron density fluctuations. Correlated variations over times long compared to the 

nominal refractive scintillation time are observed, suggesting that larger-scale density 
fluctuations are important. For these pulsars, the scintillation bandwidth as a function 
of frequency has a power-law index (~ 3.6) much less than expected for Kolmogorov 
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| turbulence (~ 4.4). Sloping fringes are commonly observed in the dynamic spectra, 

especially for PSR B0329+54. The detected range of fringe slopes are limited by our 
observing resolution. Our observations are sensitive to larger-scale fringes and hence 
, smaller refractive angles, corresponding to the central part of the scattering disk. 
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1 INTRODUCTION 

Diflractive interstellar scintillation (DISS) occurs when ra- 
diation from a pulsar has phase fluctuations imposed on it 
due to small scale variations in the interstellar electron den- 
sity. The scattered rays interfere constructively and destruc- 
tively, resulting in a modulation of the pulse intensity as a 
function of frequency and position on the observer plane. 
Pulsars are good probes of DISS because of their very small 
angular size. The relative motion of the pulsar, scattering 
material and the observer transforms the spatial variation 
into a time variation of the signal intensity. The modulation 
is characterized by a timescale (Aid) of minutes to hours and 
a decorrelation frequency scale (Ai/d) over a wide frequency 
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range from kHz to MHz, with distant pulsars scintillating 
faster in both frequency and time (Cordes, Pidwerbetsky & 
Lovelace 1986; Gupta, Rickett & Lyne 1994; Bhat, Rao & 
Gupta 1999). Longer-term fluctuations in pulsar flux densi- 
ties result from a focussing and defocussing of the ray bun- 
dle due to larger-scale inhomogeneities in the interstellar 
electron density. This phenomenon, known as refractive in- 
terstellar scintillation (RISS), is broad-band and shows less 
modulation in distant pulsars. 

By analogy with neutral gas turbulence theory, the den- 
sity fluctuations in the ionized interstellar medium (ISM) 
can be described by a power-law spectrum. The three- 
dimensional spatial power-law spectrum is of the form: 

Pan = Cl q-P (1) 

where is a measurement of the mean turbulence of elec- 
tron density along the line of sight and q = 2n/s is the 
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wavenumber associated with the spatial scale of turbulence 
s. The spatial scale s is in the range of Si nn <C s <§C s out , 
where Si nn and s ou t correspond to the inner and outer scales 
of the density fluctuations. The index /3 is thought to lie 
in the range 3 < /3 < 5. For propagation of turbulent en- 
ergy from large scales to small scales, the well-known Kol- 
mogorov theory of turbulence gives a value of /3 = 11/3. 
Observations show that this generally applies in the inter- 
stellar medium between scales of s out ~ 10 17 — 10 19 cm 
and Si nn ~ 10 8 - 10 10 cm (Spangler & Gwinn 1990; Kaspi & 
Stinebring 1992; Armstrong, Rickett & Spangler 1995; Bhat, 
Gupta & Rao 1999). 

There are discrepancies in the observed scintillation 
spectra compared to the predictions of Kolmogorov turbu- 
lence (Gupta, Rickett & Coles 1993). For example, mod- 
ulation depths and time scales often do not scale accord- 
ing to the Kolmogorov model. Various modifications to the 
Kolmogorov model have been proposed, including the steep 
spectrum model (j3 > 4) (Blandford & Narayan 1985; Good- 
man & Narayan 1985) and the inner scale model are notable. 
The steep spectrum model does not require a turbulent cas- 
cade; it could result from the superposition of noninteract- 
ing discrete structures in space. The inner scale model has a 
cutoff at the scale at which the turbulent energy dissipates 
(Coles et al. 1987). 

The variation in observed pulsar flux density as a func- 
tion of time and frequency is known as the dynamic spec- 
trum and the interference maxima in the dynamic spectrum 
are called scintles. Decorrelation time scale and bandwidth 
are parameters which quantify the scintle size. A simplified 
model is often used in which a thin scattering disk is as- 
sumed to lie midway between the pulsar and observer. In 
general this is in good agreement with observations. For a 
pulsar at distance D from the Earth, the decorrelation time 
scale for DISS has the form of 
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where v is the radio frequency (Rickett 1977; Goodman & 
Narayan 1985). In particular, for a Kolmogorov spectrum, 

AUxS'^'D-^-'V- 1 . (4) 
The decorrelation bandwidth for DISS is 

A^ d OC V W^-V D -/3/03-2) C g <4)j (5) 

Av A oc v^-V D-e^ 6 ^ (/3 > 4) (6) 

(Rickett 1977; Goodman & Narayan 1985). In the Kol- 
mogorov case, 

Avtoc^D- 11 ' 5 . (7) 

Equations 2 to 7 indicate that fluctuations are more rapid 
and have narrower bandwidths for more distant pulsars and 
at lower frequencies. 

In the Kolmogorov case, the scaling factor C\ is given 

by 



C« w 0.002 v^D-^Av- 5 ' 6 , 
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Quasi-periodic fringes resulting from interference be- 
tween multiple images of the pulsar by the scattering screen, 
are often observed in dynamic spectra (Hewish, Wolszczan 
& Graham 1985; Wolszczan & Cordes 1987; Bhat, Rao & 
Gupta 1999), with broad and narrow fringes indicating small 
and large separations of the ray paths (or images) respec- 
tively. The fringes are generally sloped on the dynamic spec- 
tra, i.e., periodic in both frequency and time and crossed 
fringes (i.e., slopes of both signs) are common. Such sloped 
fringes result from motion of the fringe pattern across the 
observer plane due to relative motion of the pulsar, Earth 
and the scattering screen. In the case of an equivalent thin 
screen located approximately midway between the pulsar 
and the observer, the maximum fringe slope is given by 

M_D6 L 

dv V e V 1 ' 

where 9 T is the refraction angle and V s is the velocity of 
the fringe pattern across the observer plane, normally dom- 
inated by the pulsar velocity (Hewish 1980; Bhat, Gupta 
& Rao 1999). The refractive angle 9 r is approximately pro- 
portional to v~ 2 for a given gradient in refractive index. 
Therefore, for large-scale gradients, the fringe slope dt/dv 
might be expected to vary approximately as v~ A . For dif- 
ferent projections of the velocity on the line between the 
images, the slope varies in magnitude and sign, resulting 
in crossed fringes of varying slope. Earlier papers discussed 
Equation 9 in terms of the frequency dependence of the re- 
fraction angle (e.g. Gupta, Rickett & Lyne, 1994), but this 
equation results directly from consideration of the geometric 
delays. For interference between a strong central image and 
other points within the midplaced scattering disk, the fringe 
rates in frequency (/„) and time (/ t ) are related by 

D 
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where v is in GHz, D is in kpc, Afd is in MHz and Cn is in 
units of m- 20/3 (Cordes 1986). 



where c is speed of light (Stinebring et al. 2001), and V B 
is dominated by the pulsar proper motion. This shows that 
crossed fringe patterns transform to parabolic arc structures 
in the 'secondary spectrum', that is, the two-dimensional 
Fourier transform of the dynamic spectrum. Such parabolic 
arcs and other more complex structures have been observed 
in a number of pulsars (Hill et al. 2003). 

In this paper, we describe observations of five pulsars 
(PSRs B0329+54, B0823+26, B1929+10, B2020+28 and 
B2021+51). In Section 2 we introduce the observations and 
data analysis. Section 3 presents the dynamic spectra and 
in Section 4 & 5 we show the auto-correlation functions and 
describe the time variations of the scintillation parameters. 
Secondary spectra for PSR B0329+54 are presented in Sec- 
tion 6. In Section 7 we discuss our results and compare them 
with previous observations and the predictions of scintilla- 
tion theories, and Section 8 summarises the results. 



2 OBSERVATIONS AND DATA ANALYSIS 

The observations were made by using the Nanshan 25-m 
telescope, operated by Urumqi Observatory, National Astro- 
nomical Observatories of China. The pulsars were regularly 
monitored from 2001 January to 2002 June with an average 
interval between observations of nine days. Our observations 
were centred at 1540 MHz, a relatively high frequency at 
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which pulsar scintillation properties have not been well stud- 
ied. Most earlier observations (e.g., Cordes, Wesberg & Bori- 
akoff 1985; Gupta, Rickett & Lyne 1994; Bhat, Rao & Gupta 
1999) were less frequent and centred on lower frequencies, 
327, 408 and 610 MHz. The receiver was a dual-polarisation 
room-temperature system with noise temperature of 95 K, 
equivalent to a system flux density of 1080 Jy. Frequency 
resolution was provided by a filterbank system consisting 
of 128 channels with channel bandwidths of 2.5 MHz giv- 
ing a total bandwidth of 320 MHz for each polarisation. 
This system was described in more detail by Wang et al. 
(2001). The data were folded online and saved to disk ev- 
ery four minutes. Total integration time was from three to 
six hours, depending on the pulsar scintillation time scale 
Aid. We made 6-hour observations for the closer pulsars, 
PSRs B1919+10, B2020+28 and B2021+51, while the obser- 
vations were three hours for the more distant pulsars PSRs 
B0329+54 and B0823+26. Between 38 and 65 observations 
were made of each pulsar. We calibrated the pulsar flux den- 
sity scale using ten strong, relatively distant pulsars which 
all have well-measured flux densities in the pulsar catalogue. 

Parameters for the observed pulsars are listed in Ta- 
ble 1. Column 3 gives the pulsar dispersion measure and 
column 4 gives the pulsar distances obtained from paral- 
lax observations. Uncertainties in the last quoted digit are 
given in parentheses. Column 5 gives the transverse veloc- 
ities based on measured proper motion and the parallax 
distance. Columns 6 and 7 give the pulsar Galactic longi- 
tude and latitude respectively, and pulsar distances from 
the Galactic plane are given in last column. 

To get the dynamic spectrum, we first dedispersed and 
summed the data for each observation to determine the pulse 
phase at the centre frequency. Pulse intensities were then de- 
termined for each sub-integration and each frequency chan- 
nel by summing the data within phase ±0.05 of the predicted 
phase of the pulse peak and subtracting a baseline offset. 

Persistent narrow-band interference was found in fre- 
quency ranges 1465-1472 MHz, 1482-1490 MHz, 1522- 
1527 MHz, 1545-1557 MHz and 1620-1625 MHz; in total 
about 15% of the frequency channels were affected. We used 
a linear interpolation across affected channels to remove the 
interference from the dynamic spectra. Occasionally, impul- 
sive broad-band interference which lasted a few minutes was 
observed; a similar interpolation was used to remove it. 



3 DYNAMIC SPECTRA 

In this section we present dynamic spectra 
for PSRs B0329+54, B0823+26, B1929+10, B2020+28 and 
B2020+51. Spectra selected to show the range of observed 
features for each pulsar are given in Fig. 1 to Fig. 5. The 
observing start time, pulsar name and MJD are given at the 
top of the each spectrum. 



3.1 PSR B0329+54 

PSR B0329+54 is the strongest northern pulsar with a dis- 
tance of 1.06 kpc. At 1540 MHz, the pulsar scintillates with 
typical time scales of 10 - 30 min and decorrelation band- 
widths of 5 - 15 MHz. We observed PSR B0329+54 at 64 



epochs over the 18-month interval, with a typical observa- 
tion time of three hours. Fig. 1 shows 12 of the dynamic 
spectra which have been corrected for both the short broad- 
band interference and the longer-term narrow-band interfer- 
ence. 

The dynamic spectra shown in Fig. 1 reveal very differ- 
ent properties compared to the lower frequency observations 
reported by Gupta, Rickett & Lyne (1994) and Stinebring, 
Faison & McKinnon (1996) which were at 408 MHz and 
610 MHz respectively. There is a much wider variation in 
scintillation time-scales and bandwidths from day to day in 
the 1540 MHz data and many spectra show multiple fringes, 
which are usually sloping, for example, the dynamic spec- 
tra obtained on MJDs 52068.6, 52078.3, 52082.1, 52131.1 
and 52202.5. Broad and narrow fringes are often seen at the 
same time, and these often have opposite slopes. Examples 
are spectra for MJDs 52027.4, 52078.3 and 52082.1. The dy- 
namic spectrum for MJD 51972.4 is interesting in that it has 
narrow vertical fringes. 

3.2 PSR B0823+26 

In total we obtained 48 dynamic spectra for PSR B0823+26, 
and Fig. 2 shows 4 of them. This pulsar has a distance of 
0.38 kpc and a large transverse velocity of 194 km s _1 . These 
parameters result in a wide decorrelation bandwidth and a 
short decorrelation time scale. As will be discussed in Sec- 
tion 7, the average scintillation scale of PSR B0823+26 is 
about 13 min and the decorrelation bandwidth is about half 
the receiver bandwidth. As seen with other pulsars the scin- 
tle time-scale is shorter than the interval between scintles at 
a given frequency. 

Despite the frequency difference between our observa- 
tions and those of Gupta, Rickett & Lyne (1994) at 408 MHz 
and Bhat, Rao & Gupta (1999) at 327 MHz, the dynamic 
spectra shown in Fig. 2 are similar with large fractional 
decorrelation bandwidths and slow drifting. The character 
of the dynamic spectra often change dramatically over short 
time intervals, for example, from a very broad single scintlc 
to multiple sloping scintles in observations on MJDs 52068.5 
and 52078.4 respectively. Clear reversals in the direction of 
drift slope are seen in observations made just one week apart, 
e.g., those on MJDs 52249.0 and 52256.0. 

3.3 PSR B1929+10 

We obtained dynamic spectra for PSR B1929+10 at 38 
epochs. This pulsar, with distance of 0.33 kpc, scintillates 
slowly and its decorrelation bandwidth is obviously compa- 
rable to or wider than the receiver bandwidth. As shown in 
Fig. 3, the six-hour observations typically were able to cover 
only one or two scintles in the time domain. Features of the 
spectra at 1540 MHz are very similar to those observed at 
lower frequencies by Gupta, Rickett & Lyne (1994). 

3.4 PSR B2020+28 

We observed 52 dynamic spectra for PSR B2020+28. At ear- 
lier stages most of the observations were about 3 hours but 
later we commenced 6-hour observations in order to cover 
several scintles in each observation. In Fig. 4 we present 4 
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Table 1. Parameters for the observed pulsars. 
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dynamic spectra. Despite its DM of 24.6 pc cm~ 3 and a dis- 
tance of 2.7 kpc, PSR B2020+28 shows rather broad scintles 
in both time and frequency. 

3.5 PSR B2021+51 

As for PSR B2020+28, PSR B2021+51 initially was ob- 
served for three hours each day; from mid-November 2001 
observation times were increased to six hours. Altogether we 
observed 65 dynamic spectra for this pulsar, and 8 examples 
are shown in Fig. 5. Clear drifting bands were observed on 
several occasions, for example, the first observed dynamic 
spectrum at MJD 52016.2. This pulsar shows two types of 
dynamic spectrum. In some cases the scintles are broad in 
both frequency and time (e.g., at MJD 52302.3), whereas in 
other cases the scintles are narrower (e.g., at MJD 52269.5). 
There appears to be a correlation between signal strength 
and scintillation type, with the pulsar generally being weaker 
when the scintles are narrow. 



4 SCINTILLATION PARAMETERS AND 
FLUX DENSITIES 

4.1 Auto-correlation functions 

The parameters of DISS can be obtained by forming a two- 
dimensional auto-correlation function (ACF) of the dynamic 
spectra S(u,t), i.e., the pulsar flux density at frequency v 
and time t. The ACF is defined by 

F(Av, At) = 

[N(Av, At)]- 1 AS ("> *) A5 (" + A ^ + At )> ( U ) 

V t 

where AS{v, t) — S{v, t) — S, and S is the mean pulsar flux 
density over each observation. N(Ais, At) is the number of 
correlated pixel pairs. The normalised ACF is then given by: 

p{Av,At) = F(Az/,Ai)/F(0,0). (12) 

We note that in Cordes (1986) the flux density 
was normalised by subtracting the mean over each sub- 
integration to minimise the effect of long-term variations, 
i.e., AS(v,t) — S(i/,t) — S(t). However this method can 
only be applied to spectra that contain many scintles, since 
in this case the sub-integration mean is close to the global 



mean apart from long-term variations. Because of our higher 
observation frequency, the number of scintles is often small 
and subtracting the mean from each sub-integration would 
affect the diffractive parameters. 

Four examples of normalised two-dimensional ACFs 
and one- dimensional cuts at zero lag are plotted in Fig. 6 
for each of the five pulsars. Following convention, the scin- 
tillation time scale Ai d is defined as the time lag at zero 
frequency lag at the point where the ACF is 1 /e of the max- 
imum, and the decorrelation frequency scale Afd is defined 
as the half- width at half-maximum of the ACF along the fre- 
quency lag axis at zero time lag (Cordes 1986). The ACFs 
show spikes at zero time and frequency lag. These spikes 
are due to noise in the dynamic spectra, and are more sig- 
nificant for weaker pulsars. Slightly wider spikes resulting 
from the interference excision process are often also present 
around zero frequency lag. To remove these spikes we used 
the neighbouring frequency-lag points to fit for a parabola 
across zero frequency lag. The central points are interpo- 
lated from this fit. 

4.2 Scintillation parameters 

We estimated the diffractive scintillation parameters for 
each observation by fitting a two-dimensional elliptical 
Gaussian function (Bhat, Rao & Gupta 1999) to the 
smoothed ACF using the non-linear least-squares fitting rou- 
tine LMM. The Gaussian function has the form: 

C g (v, t) = C exp(Ci v 2 + C 2 vt + C A t 2 ), (13) 



where Co is held fixed to unity as the ACF is normalised. 
The scintillation parameters Aid and Ava are then given by: 




Fig. 6 shows that the Gaussian function represents the 
observed ACFs well in most cases, although it may not 
be very reliable when the decorrelation bandwidth is wider 
than the receiver bandwidth, as in most observations of PSR 
B1929+10. Time variations of Aid and Afd from all obser- 
vations of the five pulsars are shown in Fig. 7 to Fig. 11 and 
average values are given in columns 2 and 3 of Table 2. One- 
sigma error estimates from the least-squares fit are plotted 
on each point in the figure. Errors in the table are one- 
sigma uncertainties in the mean values, assuming a normal 
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Figure 1. Dynamic spectra for PSR B0329+54 at center frequency 1540 MHz, with frequency and time resolution of 2.5 MHz and 4 min 
respectively. Intensity is represented by a linear greyscale between 1% (white) and 85% (black) of the maximum. 
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Figure 2. As in Fig. 1, dynamic spectra for PSR B0823+26. 



distribution. There is an additional statistical error related 
to the finite number of scintles in the dynamic spectrum. 
Following Cordes, Weisberg & Boriakoff (1985), we take the 
fractional uncertainties in At d and Af d to be given by N~ 0,5 
where N = T b s Ai/ D b s /(4 At d Av d ) is the number of observed 
scintles. Typical values of these additional uncertainties are 



shown by error bars in the upper left corner of the Atd and 
Ais d plots in Fig. 7 to 11. 

Fitting for the two-dimensional Gaussian ellipse pro- 
vides parameters for the sloping features of the dynamic 
spectrum. Following Bhat, Rao & Gupta (1999), we adopt 
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Figure 3. As in Fig. 1, dynamic spectra for PSR B1929+10. Note the different time scale on the vertical axis compared to Fig. 1. 
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Figure 4. As in Fig. 1, dynamic spectra for PSR B2020+28. 



dt/dv as the drift rate; in terms of the fitted parameters 
dt/dv is given by: 

t-m- 

The slope visibility r is given by: 

r = C 2/ V4Ci C 3 . (16) 

Derived values for the drift rate and absolute value of slope 
visibility (|r|) are plotted in Fig. 7 to Fig. 11 and average 
values are given in columns 5 and 6 of Table 2. 

The scattering strength u is defined to be the ratio of 
Fresnel scale sf to the DISS scale Sd (Rickett 1990; Gupta 
1995); u ^> 1 corresponds to strong scattering and u <§C 1 
to weak scattering. In terms of observable quantities, u is 
expressed as 



(17) 



(Bhat, Gupta & Rao 1999), where v is the observing fre- 
quency. Observed mean values of u are given in column 7 of 
Table 2. Scattering is stronger for PSR B0329+54 compared 
to PSRs B0823+26, B2020+28, and B2021+51, despite their 
similar DMs. 

Based on the simple thin screen model with power-law 
density fluctuations, the timescale for refractive scintillation 
t r is given by 



2v 
Av~ d 



(18) 



(Rickett 1990). This parameter is given in column 8 of Ta- 
ble 2. We use Equation 8 to estimate mean values of the 
fluctuation scaling parameter C„; derived values are given 
in column 9 of Table 2. 

Pulsar velocities can be estimated from the scintillation 
parameters Aid and Ai^d. For a thin scattering screen placed 
midway between the pulsar and the observer, the transverse 
velocity of the pattern is given by: 



V a « 3.85 x 10 



4 VAv d D 
vAt d 



(19) 



where V B is in km s _1 , A^d is in MHz, the pulsar distance D 
is in kpc, v is in GHz, and Aid is in seconds (Gupta, Rickett 
& Lyne 1994). V a is a combination of the transverse velocity 
of the pulsar, the scattering medium and the Earth's mo- 
tion. The latter two components are usually comparatively 
small, so V s represents the transverse velocity of the pulsar. 
From Equation 19 and the parallax distances given in Ta- 
ble 1, we derived transverse velocities for the five pulsars. 
The resulting values are plotted in Fig. 7 to Fig. 11 and 
average values are given in column 10 of Table 2. 
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Figure 5. As in Fig. 1, dynamic spectra for PSR B2021+51. 



4.3 Flux densities 

Long-term observations have shown that pulsar flux density 
(S) variations are dominated by refractive scintillation. Dis- 
tant pulsars show little flux density variation, implying that 
intrinsic pulsar radio luminosities are stable (Kaspi & Stine- 
bring 1992). Observed flux densities are plotted in Fig. 7 to 
Fig. 11 and mean values are given in column 11 of Table 2. 



5 MODELLING THE ACFS 

The observations of PSR B0329+54 gave us 64 independent 
estimates of the frequency-time correlation function of inten- 
sity. Even though these often show drifting features the aver- 
age correlation function should still be estimated with good 
confidence and we compared it with its theoretical shape 
(assuming strong scattering) from a thin scattering medium 
with a Kolmogorov spectrum. 

The theoretical intensity correlation versus spatial and 
frequency offsets, assuming a Kolmogorov spectrum and a 
thin screen, has been derived by Lambert & Rickett (1999; 
see their Figure 5) . In Fig. 12 we compare their model with 
the mean of the observed auto-correlation functions, by as- 
suming that time offsets translate simply into spatial offsets. 
We estimated Au^ by fitting the theoretical shape of the cut 
along the frequency axis and Atd from the cut along the time 
axis. The observed one-dimensional correlation functions are 
shown by the error bars (three times the nominal standard 
error in the mean) with a solid line for the model. 



These plots show good general agreement in shape along 
both axes. Particularly evident is the steeper slope near zero 
lag in the frequency domain than in the time domain for 
both the model and the observations. The minimum chi- 
squared value is about twice the number of data points in the 
fit, which maybe due to too low an estimate of the error in 
the mean correlation function. We suspect that the temporal 
variations in the two-dimensional correlations do not follow 
normal statistics which causes the underestimate of the error 
in its mean. 

Whereas along the axes the model agrees with the the- 
ory, away from the axes the theory contours bulge outward 
more than the observed contours. The outward bulging may 
be related to the arc phenomenon in the secondary spectrum 
(see Section 6). We also note that scattering extended along 
the line of sight causes a reduction in the bulging effect, as 
can be seen in Figure 5 of Lambert & Rickett (1999). 



6 SECONDARY SPECTRA 

Secondary spectra, that is, two-dimensional Fourier trans- 
forms of dynamic spectra, give additional insight into scin- 
tillation and scattering in the interstellar medium. Fig. 13 
shows secondary spectra, for several observations of PSR 
B0329+54 where significant fringing was observed. These 
spectra are much more amorphous and centrally concen- 
trated than those shown by Hill et al. (2003) for other pul- 
sars at a similar frequency. We note that their observations 
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Figure 6. Samples of normalised two-dimensional auto-correlation functions of the dynamic spectra shown in Fig. 1 to Fig. 5. The 
contour interval is 0.1, with dashed lines representing negative contour levels. The solid line in the two smaller plots at the bottom 
and right are the one-dimensional ACFs at zero lag in time and frequency respectively. The dotted lines are zero-lag cuts through the 
two-dimensional Gaussian fits to the observed data. 
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Table 2. Average scintillation parameters of the observed pulsars. 
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Figure 7. Time variations of At<j, Ai/j, main axial angle, 
slope visibility, transverse velocity and the flux density for PSR 
B0329+54. Error bars (±1<t) plotted on each point arc from the 
2D-clliptical Gaussian fitting. The statistical uncertainty related 
to the finite number of observed scintles is indicated by the ±lcr 
error bar in the mean value shown at the top-left corner of the 
Atd an d A^j plots. 



have narrower bandwidths and shorter observation times but 
much higher time and frequency resolution (10 - 30 s and 
~0.1 MHz respectively) compared to 4 min and 2.5 MHz for 
our observations. Our observations therefore are sensitive to 
larger-scale fringes and hence smaller values of 9 r than those 
of Hill et al. (2003), i.e., we are observing the central part 
of the scattering disk with high spatial resolution, whereas 
the higher-resolution observations are senstive to large re- 
fraction angles and very fine fringing structure. Although 
there are some indications of rudimentary arc structure in 
Fig. 13, the model of interference between a strong central 
image and surrounding discrete images is clearly less appro- 
priate — we are in fact observing the structure of the central 
image. 

Across our wide observed bandwidth (A^ obs = 
320 MHz), the frequency dependence of the slope (dt/dv oc 
v~ A ), sec Equation 9) corresponds to a 2:1 variation with 
larger slopes (more vertical fringes) at lower frequencies. 
There is some evidence for this in Fig. 1, for example, in the 
narrow fringing visible in the spectrum for MJD 52202.5. 
This slope variation will smear the corresponding feature 
in the secondary spectrum, contributing to the amorphous 
appearance of the secondary spectra (Fig. 13). 



7 DISCUSSION 

7.1 Dynamic spectra 

For PSR B0329+54, the character of the dynamic spectra 
is very different from that observed at lower frequencies. 
Sloping fringes of different frequencies and slopes are very 
common in the 1540 MHz observations, whereas at lower 
frequencies they are uncommon (e.g., Stinebring, Faison & 
McKinnon, 1996). The frequency dependence of decorrela- 
tion bandwidth is directly observable within our receiver 
bandwidth in many dynamic spectra. Taking the observa- 
tion at MJD 52294.5 for PSR B0329+54 as an example, the 
width of the scintle at 1400 MHz is about 13 MHz, whereas 
at 1640 MHz it is about 27 MHz, consistent with the pre- 
diction from Equation 7. 

Despite PSRs B0329+54, B0823+26, B2020+28 and 
B2021+51 having very similar dispersion measures, their 
scintillation properties are quite different. As is well known, 
Galactic latitude appears to be important with stronger 
scattering for PSR B0329+54. Distance seems less impor- 
tant with very different DISS properties for PSR B0823+26 
and PSR B1929+10. 
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Figure 8. As in Fig. 7, time variations of Aid, A;/d, main axial 
angle, slope visibility, transverse velocity and the flux density for 
PSR B0823+26. 
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Figure 9. As in Fig. 7, time variations of Aid, Ai/d, main axial 
angle, slope visibility, transverse velocity and the flux density for 
PSR B1929+10. 



7.2 Time variations 

Observed dynamic spectra vary greatly from observation to 
observation, in all respects: decorrelation time scale, band- 
width, brightness and the characteristics of sloping features. 
Fig. 7 to Fig. 11 show that scintillation parameters have 
a wide range of variation. For example, for PSR B0329+54 
(Fig. 7), the observed diffractive timescale and decorrelation 
bandwidth vary over the ranges 12-30 min and 5-34 MHz 
respectively. Similar variations are observed for the other 
four pulsars. The interval between observations is typically 
several times the refractive timescale t r (Table 2) and so we 
would expect uncorrelated variations. However, the expected 
modulation due to refractive scintillation can only account 
for a small amount of the observed variations. Modulation 
indices (ratio of rms deviation to mean value) for the varia- 
tions in Atd, Az/ d and flux density, m t , m b and m r respec- 
tively, are given in Table 3. Also given in the Table are pre- 



dicted modulation indices for the case of a thin screen with 
a Kolmogorov density spectrum (f3 = 11/3) and strong scat- 
tering (Romani, Narayan & Blandford 1986; Bhat, Gupta & 
Rao 1999). In all cases, the observed modulation indices are 
greater than the predicted values. With our relatively high 
observing frequency and relatively nearby pulsars, the scat- 
tering strength u is only modestly greater than one, so the 
strong scattering assumption is not really satisfied. Also, the 
form of the variations suggests that large-scale fluctuations 
in the screen are important, so the fluctuation spectrum is 
likely to be steeper than in the Kolmogorov case. 

Fig. 7 to Fig. 11 show that there are significant cor- 
related changes in decorrelation time-scale, bandwidth and 
drift rate over intervals of several weeks and in some cases 
many months, much larger than the predicted refractive 
timescale £ r . For example, in Fig. 7 for PSR B0329+54, there 
are high values of decorrelation bandwidth at MJD 52035 
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Figure 10. As in Fig. 7, time variations of Atd, Ai/j, main axial 
angle, slope visibility, transverse velocity and the flux density for 
PSR B2020+28. 
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Figure 11. As in Fig. 7, time variations of Atd, A^d, main axial 
angle, slope visibility, transverse velocity and the flux density for 
PSR B2021+51. 



and 52287, low values of drift rate between MJD 52295 and 
52370, an increase and then decrease in slope visibility be- 
tween MJD 52160 and 52215 and larger than average flux 
densities between 52050 and 52110. These results show that 
gradients in refractive index of scale much larger than that of 
the diffractive scattering disk exist in the scattering screen. 

In general, there is little correlation between variations 
in the various parameters. However, there are exceptions. 
For example, Fig. 11 shows a very significant correlation be- 
tween decorrelation time-scale and flux density. As remarked 
in Section 4.2, the pulsar tends to be strong when the scin- 
tillation timescale is long. This suggests that most of energy 
from the pulsar is confined to the central region of the scat- 
tering disk where the scattering angles are small, leading to 
broad scintles. Only when we are in a null of these broad 
patterns do we see weaker radiation from outer parts of the 



scattering disk where scattering angles are larger and scin- 
tles are smaller. 

Due to the wide range variation of diffractive timescale 
and decorrelation bandwidth, the derived velocity is also 
highly variable. However, there is remarkably good agree- 
ment of the average scintillation velocities (Table 2) with the 
velocities derived from proper motion measurements (Ta- 
ble 1). We attribute this to our long-term observations which 
averaged over the refractive scintillation modulations (Bhat, 
Rao & Gupta 1999). This agreement also indicates that lo- 
cating a thin screen approximately half-way between the 
pulsar and observer is a good assumption (Gupta, Rickett 
& Lyne 1994; Nicastro et al. 2001). 
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Figure 12. Auto correlation versus frequency and time lag for 
PSR B0329+54. Solid contours for the average auto correlation of 
64 individual observations. Dashed lines are contours for a Kol- 
mogorov model in the diffractive limit of scattering from a thin 
screen. Av& and At<j were fitted to the cuts along the axes as dis- 
played in the two small bottom and right plots, where the data are 
represented by error bars that are three times the nominal stan- 
dard error in the mean. Notice how the dashed contours agree 
with the solid contours on each axis, but bulge out away from the 
axes. 



Table 3. Modulation indices for At^, Av^ and flux density. 
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7.3 Power spectrum of electron density 
turbulence 

The frequency dependence of the DISS parameters can be 
used to estimate the spectral slope of interstellar turbulence. 
Table 4 shows the results from earlier low-frequency observa- 
tions, with the data in columns 2 and 3 derived from obser- 
vations using the Ooty Radio Telescope at 327 MHz by Bhat, 
Rao & Gupta (1999), columns 4 and 5 are from 408 MHz 
observations by Gupta, Rickett & Lyne (1994) and columns 
6 and 7 are from 610 MHz observations by Stinebring, Fai- 
son & McKinnon (1996). For convenience, columns 8 and 9 
are from our observations, repeated from Table 2. For Kol- 
mogorov turbulence, the frequency dependences of scintil- 
lation parameters are given by: Atd oc v 1 ' 2 , Av^ oc u 4A 
(see Equations 4 and 7). Scintillation timescales are in rea- 
sonable accord with the predicted frequency dependence, 
but the observed values of the decorrelation bandwidth are 
systematically lower at high frequencies. Fitting for the ob- 
served time scale and decorrelation bandwidth we obtained 
the power index, a t and a v and the corresponding density 
spectral indices f3 t and p v , using a model for a density fluctu- 
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Figure 13. Secondary spectra for selected observations of PSR 
B0329+54. The spectra have been Hanning weighted to reduce 
sidelobes and improve the signal/noise ratio. The dotted lines 
indicate expected arcs derived from Equation 8. 



ation spectrum with f3 < 4 (Equations 2 & 5). These values, 
shown in the last four columns in Table 4 (except in one 
case) show a flatter frequency dependence for Atd and Avd 
than predicted, implying a density fluctuation index much 
larger than the Kolmogorov value and in fact greater than 
the critical value of four. We therefore applied a model for a 
steeper density fluctuation spectrum, j3 > 4 (Equations 3 & 
6). However, the derived values contradict the assumption 
giving /3t < 4 and /3„ < 4. 

The flat frequency dependence of At& and Av& is not 
predicted by either theory. Similar low values of a„ were also 
observed earlier by Johnston, Nicastro & Koribalski (1998) 
and Lohmer et al. (2001). The former work consists of multi- 
frequency DISS observations for large sample of pulsars and 
in the latter paper the index was derived from scatter broad- 
ening observations. This result may imply a relatively small 
outer scale for the scattering fluctuations or may be a con- 
sequence of weaker scattering at the higher frequencies. 

By comparing the diffractive and refractive scattering 
angles, 8d and 6 V respectively, we can estimate the spectral 
index /3 of the density fluctuations (Equation 1). In terms 
of observable quantities, they are given by 
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Table 4. Fitting for frequency dependences of Atd and Ai^ and the corresponding density spectra. 
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a Bhat, Rao & Gupta 1999 
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c Stincbring, Faison & McKinnon 1996 

d Our Results 



and Equation 9 rearranged 



% = (^) fa 



\du) 



(21) 



where V s is given in Table 2 and dt/du is the drift slope of 
the dynamic spectrum from Equation 15. The mean values 
of (6 d ) and (6i) are given in columns 2 and 3 of Table 5. As 
expected for random fluctuations in the scattering medium, 
(6 T ) is essentially zero within the uncertainties. However, the 
scale of the refractive scattering can be estimated from the 
rms fluctuation of 8 T , denoted as og T and given in column 4 
of Table 5. Following Bhat, Gupta & Rao (1999) we scale 
this quantity by a factor of v2 to allow for mis-alignment 
of the refractive gradient and the pattern velocity. 

We estimate the spectral index of the electron density 
fluctuations using 



/3 = 4 + 



log(a 9r /<0 d » 



logtt 



(22) 



where u is the measurement of scattering strength defined 
in Equation 17 (Bhat, Gupta & Rao 1999). This formula 
indicates that ae I < (8 d ) corresponds to a flatter spectrum 
f3 < 4 but ae, > (9 d ) corresponds to a steeper spectrum 
P > 4. In our case as a Bl <C (# d ), and the measured /3 are 
in the range 3.40 - 3.75, on average somewhat less than 
the Kolmogorov spectrum index of 3.67, but in most cases 
consistent with it. 

Taken at face value, these results are inconsistent with 
those based on the frequency dependence of scintillation pa- 
rameters (Table 4). One possible explanation is that the 
fluctuations at the diffractive and refractive scales are in- 
dependent, not related by the cascade of turbulent energy 
and hence not described by the Kolmogorov spectrum (Arm- 
strong, Rickett & Spangler 1995). We also note that our 
range of derived drift slopes is limited by the observing pa- 
rameters, especially time and frequency resolution. Higher 
resolution observations (e.g. Hill et al. 2003) detect a wider 
range of slopes and hence imply larger refractive angles and 
a steeper spectrum for the density fluctuations. 



Table 5. Estimates of diffractive and refractive 
scattering angles, and the spectral index of elec- 
tron density fluctuation. 
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3.59 


0823+26 
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0.004(5) 


0.044 


3.75 


1929+10 


0.043(2) 


0.003(2) 


0.019 


3.33 


2020+28 


0.029(1) 


-0.0003(9) 


0.009 


3.40 


2021+54 


0.039(1) 


-0.002(2) 


0.017 


3.61 



8 SUMMARY 

In this paper we present dynamic spec- 
tra for PSRs B0329+54, B0823+26, B1929+10, B2020+28 
and B2021+51 over a 320 MHz centred at 1540 MHz. The 
time and frequency bandwidth resolutions are 4 min and 
2.5 MHz respectively. The dynamic spectra show different 
features compared to earlier lower- frequency observations. 
For example, PSR B0329+54 presents frequent drifting pat- 
terns and fringes which were not seen at lower frequencies. 
Our individual observations are generally long enough to 
cover at least a few scintles, but in some cases, especially 
for PSR B1929+10, the scintles are wide compared to our 
bandwidth. In these cases, the value of Afd is based on just 
the peak of the ACF and is rather uncertain. Although the 
observations show a wide range of variations in diffractive 
scintillation parameters, the average scintillation velocities 
are in good agreement with the proper motion velocities. 

The frequency and time correlations are successfully de- 
scribed by theory for a thin Kolmogorov scattering screen. 
However, there is a consistent disagreement for points off the 
axes, which suggests turbulence extended along the line of 
sight. There is also much greater variability than expected 
in the shapes of the daily auto-correlations, which suggests 
a stronger than expected influence of scales larger than the 
scattering disc. We note here that similar conclusion has 
been drawn by Lambert & Rickett (2000), Spangler (2001) 
and Shishov et al. (2003). These authors also pointed to 
more complex models than a Kolmogorov thin screen. We 
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do not have theoretical predictions for the full frequency- 
time correlation for the implied steeper spectrum models, 
which thus becomes an area for future work. 

With our wide observed bandwith, long observation 
times and lower time and frequency resolution, our observa- 
tions are sensitive to larger-scale fringes and hence smaller 
refractive angles, corresponding to the central part of the 
scattering disk. The secondary spectra we obtained for PSR 
B0329+54 are centrally concentrated and rather amorphous. 
Some of the smearing results from a v~ i frequency depen- 
dence of fringe slope across our wide observed band. 

Modulations in time scale, decorrelation bandwidth and 
flux density are greater than predictions from a thin screen 
model with a Kolmogorov spectrum of density fluctua- 
tions. The long-term variations in DISS parameters sug- 
gests larger-scale density fluctuations exist, corresponding 
to a steeper fluctuation spectrum. Comparison of the ob- 
served DISS timescales with values extrapolated from lower 
frequency measurements based on the Kolmogorov spec- 
trum show good agreement. However, large discrepancies 
were found in extrapolated decorrelation band widths, with 
the observed values substantially smaller than the predicted 
values. This gives a flatter frequency dependence of decor- 
relation bandwidth, corresponding to a steeper density fluc- 
tuation spectrum. 

In contrast, estimation of the power-law density fluctua- 
tion index /3 from diffractive and refractive scattering angles 
gives a value in the range 3.40 - 3.75, on average smaller but 
consistent with a Kolmogorov spectrum. The range of de- 
rived drift slopes is limited by the observing parameters, and 
observations with higher resolution will detect larger refrac- 
tive angles and hence imply a steeper density spectrum. It is 
also possible that the large-scale fluctuations responsible for 
refractive effects are not directly related to the smaller-scale 
fluctuations responsible for diffractive scintillation. 
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